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— , ■ Abstract 

The stationary, spherically symmetric accretion of dark energy onto a Schwarzschild 
black hole is considered in terms of relativistic hydrodynamics. The approximation of 
an ideal fluid is used to model the dark energy. General expressions are derived for the 
qq 1 accretion rate of an ideal fluid with an arbitrary equation of state p = p(p) onto a black 

I ' hole. The black hole mass was found to decrease for the accretion of phantom energy. The 

\Q , accretion process is studied in detail for two dark energy models that admit an analytical 

' solution: a model with a linear equation of state, p — a(p — po), and a Chaplygin gas. For 

one of the special cases of a linear equation of state, an analytical expression is derived for 
the accretion rate of dark energy onto a moving and rotating black hole. The masses of 
all black holes are shown to approach zero in cosmological models with phantom energy in 
pH , which the Big Rip scenario is realized. 

1 INTRODUCTION 

In recent years, strong observational evidence that the Universe is currently expanding with 
acceleration has been obtained. In the Einstein theory of gravitation, this positive acceleration 
is explained by the dominance of dark energy with a negative pressure in the Universe [1- 
4]. Several theoretical models of dark energy have been suggested: the vacuum energy (the 
cosmological constant A) or such dynamical components as quintessence [5-10] and k-essence 
[11-13], Models with dynamical dark energy seem more realistic, since tracker [14, 15], or 
attractor, solutions are realized in them. Thus, the problem of fine tuning the parameters of 
the Universe is solved [11-13], 

A peculiar property of cosmological models with dark energy is the possibility of a Big 
Rip [16, 17]: an infinite increase in the scale factor of the Universe in a finite time. The Big 
Rip scenario is realized in the case of dark energy, the so-called phantom energy (for which 
p + p < 0). In the Big Rip scenario, the cosmological phantom energy density tends to infinity, 
and all of the bound objects are torn apart up to subnuclear scales. It should be noted, however, 
that the condition p + p < alone is not enough for the Big Rip scenario to be realized [18] . 
Alam et al. [19] analyzed data on distant supernovas in a model-independent way and showed 
that the presence of phantom energy with —1.2 < w < —1 in the Universe at present is highly 
likely. The quantum properties of the phantom energy in curved spacetime were considered in 
[20]. The entropy of the Universe filled with phantom energy was discussed in [21]. Models with 
phantom energy are also used to construct mole burrows [22, 23]. The accretion of a scalar field 
onto a black hole from special poten-tials V{4>) was considered in [24-29]. We use a different 
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approach to describe the accretion of dark energy onto a black hole; more specifically, we model 
the dark energy by an ideal fluid with a negative pressure. 

In our recent paper [30] (see also [31]), we showed that the masses of all black holes in the 
Universe with phantom energy gradually decrease, and the black holes disappear completely by 
the Big Rip. In this paper, we consider in detail the stationary spherical accretion of dynamical 
dark energy onto a black hole. The dark energy is modeled by an ideal fluid with a negative 
pressure. The history of research on the accre-tion of an ideal fluid onto a compact object begins 
with Bondi's classic paper [32]. A relativistic generalization was made by Michel [33] (see also 
[34-41] for further generalizations and supplements to Michel's solution). Carr and Hawking [42] 
considered the accretion of dust and radiation onto a black hole by solving the complete system 
of Einstein equations and taking into account the back reaction of the surrounding matter (see 
also [43] for a description of the progress made in this area and for a discussion of fundamental 
questions). Below, we obtain a solution for the stationary accretion of a test relativistic ideal 
fluid with an arbitrary equation of state p(p) onto a Schwarzschild black hole. Using this solu- 
tion, we show that the black hole mass decreases during the accretion of phantom energy. The 
masses of black holes can decrease during accretion in the case of phan-tom energy due to the 
violation of the energy domi-nance condition (p + p > 0) that underlies the theorem on the 
nondecreasing area of the event horizon of a classical black hole [44] . 

This paper is structured as follows. In Section 2, we derive general equations for the spherical 
accretion of an ideal fluid and describe basic parameters of the steady energy flux onto a black 
hole. We consider an arbitrary equation of state, w = p/p, where the pressure p can be positive 
(for ordinary matter) and negative (for dark energy, including phantom energy w < —1). Note 
that the parameter w of the equation of state need not be constant in our approach. Accretion 
causes the black hole mass to change: the mass increases for p + p > and decreases for 
p + p < 0. The energy flux turns out to be completely determined by the black hole mass M, 
the dark energy density at infinity poo, and the equation of state p = p(p) only if < dp/dp < 1. 
In this case, there is a critical point that fixes the flux just as for an ordi-nary fluid. When 
the condition < dp/dp < 1 is violated, the dark energy flux onto a black hole can formally 
be arbitrary. For < dp/ dp < 1, we describe the method of calculating the fluid parameters 
at the critical point and the energy flux onto a black hole for given M, p^, and p — p(p). 
In Section 3, we consider specific models of the equation of state for dark energy In the first 
model, we use a simple equation of state with a linear density dependence of the pressure. We 
consider the special cases of accretion of several types of ideal fluid: ther-mal radiation, matter 
with an ultrahard equation of state, dark energy with dp/dp §s 0, and linear phantom energy. 
The accretion rate of dark energy onto a moving black hole was calculated for the special case 
of dp/ dp — 1. As the second model, we investigate the accretion of a Chaplygin gas onto a 
black hole. The evolution of the black hole mass in the Universe with the Big Rip is considered 
in Section 4. The possibility that the presence of phantom energy will lead the Universe to 
the Big Rip in the future has been discussed in recent years. The problem of the fate of black 
holes in this Universe is solved in a rather unexpected way: black holes are not torn apart, but 
disappear by the Big Rip due to the accretion of phantom energy, irrespective of their initial 
masses. In Section 5, we discuss the cor-respondence between the accretion of dark energy 
modeled by an ideal fluid onto a black hole and the accretion of a scalar field. The results 
obtained are briefly discussed in Section 6. 

2 GENERAL EQUATIONS FOR SPHERICAL ACCRE- 
TION 

Let us consider the stationary, spherically symmetric accretion of an ideal fluid that models the 
dark energy in the special case of a negative pressure onto a black hole. The dark energy density 
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is assumed to be low enough for the metric to be a Schwarzschild one with a high accuracy: 



_ / 2M\ „> / . 2 1/ 



ds 2 = (1 J eft 2 - I 1 1 dr 2 -r 2 {d6 2 + sin 2 Od<j) 2 ), (1) 

where M is the mass of the black hole, r is the radial coordinate, and 9 and 4> are the angular 
spherical coordinates. We model the dark energy by an ideal fluid with the energy-momentum 
tensor 

T^ u = (p + PtUpUv - pgpv, (2) 

where p is the density, p is the dark energy pressure, and = dx^/ds is the radial 4- velocity 
component. The pressure is assumed to be an arbitrary function of the density, p = p{p). Inte- 
grating the zeroth (time) compo-nent of the conservation law = yields the first integral 
of motion for stationary, spherically symmetric accretion (Bernoulli's rclativistic equation or 
the energy equation): 

/ 2 \ 1/2 
{p + p) f l--+u 2 \ x 2 u = d, (3) 

where x = r/M, u = dr/ds, and C\ is the constant determined below. To find the second 
integral of motion, we use the equation for the component of the energy-momentum tensor 
conservation law along the 4- velocity 



u 



T» v . v = 0. (4) 



In our case, this quation is [45] 

u^p !fl + (p + p)u^ = 0. (5) 

For the given equation of state 

p = p(p), (6) 
the auxiliary function n can be defined by the relation 

dp dn 



p + p 



(7) 



The function n is identical to the particle concentration for an atomic gas, but it can also be 
used to describe a continuous medium that does not consist of any particles. In this case, the 
'concentration' n is a formal auxiliary function. For an arbitrary equation of state p = p{p), we 
obtain a solution for n from Eq. (JZJ): 




v Poo 

Using (JSJ, we find the sought second integral of motion from Eq. 



-ux 



= -A, (9) 



where rioa (the dark energy 'concentration' at infinity) was introduced for convenience. In the 
case of a fluid flow directed toward the black hole, u = (dr/ds) < 0, and therefore, the numerical 
constant A > 0. From © and 10, we can easily obtain 
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1--+U 2 ) =C 2 , (10) 



3 



where 

n Poo+P(Poo) ni » 
"-(Poo) 

Let us now calculate the radial 4-velocity component and the fluid density on the event 
horizon of the black hole, r = 2M. Setting x = 2, we obtain from Eqs. @, l(TU|) and ltTT|l 

A pn + p(ph) = n 2 (p H ) 

4 Poo +P{Poo) ^(Poo)' 

where pn is the density on the x — 2 horizon. Thus, having specified the density at infinity p^, 
the equation of state p = p(p) , and the flux A and using definition J7J) of the concentration, we 
can calculate the fluid density pn on the event horizon of the black hole from l|12Jl . Given the 
density on the horizon pn, we can easily determine the radial fluid 4-velocity component on the 
horizon from 0: 

A n( Poo ) 

"H = — -r — -, r- (13) 

4 n(pn) 

Below, we will see that the constant A, which defines the energy flux onto the black hole, 
can be calculated for hydrodynamically stable ideal fluids with dp/ dp > 0. This can be done 
by determining the fluid parameters at the critical point. Following Michel [33], we find the 
relationship between the parameters at the critical point: 

*-t£s- (14) 



where 



Together with Q, this yields 



V 2 = n -1, (15) 

p + p an 



V' = ci(p), (16) 

where c 2 — dp/ dp is the square of the effective speed of sound in the medium. We derive the 
following relation for the critical point from Eqs. (|14(l . (|lf>|t . (|1 1|) . and (|1C)|) : 

P* + p(p*) n ,o2/ m 1/2 Poo±p{poo) n7 , 

which fixes the fluid density at the critical point for an arbitrary equation of state p = p(p). 
Specifying and using ©, we can determine n(p*). Accordingly, the quantities and u* can 
be calculated from (|14|l and . As a result, the numerical constant A can be calculated by 
substituting the derived quantities into ©. For c 2 s < or c 2 s > 1, no critical point exists beyond 
the event horizon of the black hole (a;* > 1), implying that the dark energy flux onto the black 
hole depends on the initial conditions for an unstable ideal fluid c 2 < or a 'superluminal' fluid 
c 2 > 1. This result has a simple physical interpretation: the accreted fluid has a critical point 
if its speed increases from subsonic to supersonic values as it approaches the black hole. In 
contrast, for <? s or c 2 > 1, the critical point either does not exist or is formally within the event 
horizon of the black hole. It should also be noted that fluids with c 2 < are hydrodynamically 
unstable (see [46, 47] for a discussion). 

Equation (|10fl . together with (j^Jl, ©, and @, defines the accretion rate onto a black hole. 
These equations are valid for an ideal fluid with an arbitrary equation of state p = p(p), in 
particular, for a gas of massless particles (thermalized radiation) and a gas of massive particles. 
For a gas of massive particles, the system of equations © and l|l(J|l reduces to a similar system 
of equations found by Michel [33]. It should be noted, however, that Eqs. JBJ, (JSJ), (J5J), and 
(|1L)|) are also valid for dark energy, including phantom energy with p + p < 0. In these cases, 
the concentration n(p) is positive for any p, while the constant C2 in Eq. (|ll)f) is negative. 
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The rate of change in the black hole mass (the energy flux onto the black hole) through 
accretion is 

M = -47rr 2 T r . 
Using (9) and (10), this expression can be rewritten as [30] 

M ^^AM 2 [ Poc +p{ Poc )]. (18) 

It follows from Eq. (18) than the mass of the black hole increases as it accretes the gas of 
particles when p > 0, but decreases as it accretes the phantom energy when p + p < 0. In 
particular, this implies that the black hole masses in the Universe filled with phantom energy 
must decrease. This result is general in nature. It does not depend on the specific form of the 
equation of state p = p(p); only the satisfaction of the condition p + p < is important. The 
physical cause of the decrease in the black hole mass is as follows: the phantom energy falls 
to the black hole, but the energy flux associated with this fall is directed away from the black 
hole. If we ignore the cosmological evolution of the density p^, then we find the law of change 
in the black hole mass from (18) to be 

M = Mi (l-^) , (19) 

where Mi is the initial mass of the black hole, and r is the evolution time scale: 

T = l/{4nAM i [p 00 +p(p 00 )]}. (20) 



3 ANALYTICAL ACCRETION MODELS 

3.1 Model of a Linear Equation of State 

Let us consider the model of dark energy with a linear density dependence of pressure [30] : 

p = a(p-p ), (21) 

where a and po are constants. Among the other cases, this model describes an ultrarelativistic 
gas {p = p/3), a gas with an ultrahard equation of state (p = p), and the simplest model of 
dark energy (p = and a < 0). The quantity a is related to the parameter w = p/p of the 
equation of state by w = a(p — Po)/p- 

An equation of state with w — const < throughout the cosmological evolution is com- 
monly used to analyze cosmological models. The matter with such an equation of state is 
hydrodynamically unstable and can exist only for a short period. Our equation of state (21) for 
a > does not have this shortcoming. For a > 0, it also allows the case of hydrodynamically 
stable phantom energy to be described, which is not possible when using an equation of state 
with w — const < — 1. In the real Universe, the equation of state changes with time (i.e., w 
depends on t). Therefore, Eq. (21) has the meaning of an approximation to the true equation 
of state only in a limited p range. From the physical point of view, the condition p > must 
be satisfied for any equation of state in a comoving frame of reference. In particular, the state 
of matter with p = 0, but p ^ 0, is physically unacceptable. The corresponding constraints for 
the equation of state (21) are specified by conditions (29) and (30) given below. 

For a < 0, there is no critical point for the accreted fluid flow. For a > 0, using (14) and 
(16), we obtain the parameters of the critical point 

1 + 3a o a 
X * = — U * = —a (22) 
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Note that the parameters of the critical point (22) in the linear model (21) are determined 
only by dp/ dp = a and do not depend on p , which fixes the physical nature of the fluid 
under consideration: a relativistic gas, dark energy, or phantom energy. Note also that no 
critical point exists beyond the event horizon of the black hole for a > 1 (this corresponds to 
a nonphysical situation with a supcrluminal speed of sound). 

Let us calculate the constant A, which defines the energy flux onto the black hole. We find 
from (8) that 

(23) 



Petf 



Peff,c 



where we introduced the effective density 



Peff = P + P = -Pa® + (1 + a)p 



Using (17), we obtain 



Peg, 
Peff.oc 



a/(l+a) 



(l + Sa) 1 / 2 , 



(24) 



where p c g; and p e ff,oo, are the effective densities at the critical point and at infinity, respectively. 
Substituting (24) into (23) and using (9), we obtain for the linear model 



(l + 3a)( 1+3 »>/ 2 " 
4a 3 / 2 ' 



(25) 



It is easy to see that A > 4 for < a < 1. A = 4 for a = 1 (this corresponds to c s — 1); i.e., 
the constant A is on the order of 1 for relativistic speeds of sound. Using (25), we obtain from 
(20) 

(l + 3a)( 1+3 «)/ 2 ° 



7rM i (p 0O +Poo)- 



,3/2 



(26) 



To determine the fluid density on the event horizon of the black hole, we substitute (23) into 
(12) to yield 

apo . ( apo \ ( A \ 
= t— + I Poo - t— t ' ( 2? ) 



1 + a 



1 



where A is given by (25). For < a < 1, the effective density on the horizon p e ff,H cannot be 
lower than p c ff oo. The radial 4-velocity component on the horizon can be found from (13) and 
(27): 



A 

"i. = - I - 



-a/(l-a) 



(28) 



The value of «h changes from 1 to 1/2 for < a < 1. 

The linear model (21) describes the phantom energy when p^,/ p n < a/(l + a). In this case, 
p + p < 0. However, the requirement that the density p be nonnegative should be taken into 
account. This parameter can formally be negative in the range < a < 1. Such a nonphysical 
situation imposes a constraint on the linear model (21) under consideration. For a physically 
proper description of the accretion process, we must require that the density p be nonnegative. 
We obtain the following constraint on the validity range of the linear model from (27) for 
hydrodynamically stable phantom energy: 



1 + a 



A 



-(l+a)/(l-a) 



Poo a 
po 1 + a' 



As follows from (29), at a given a, we can always choose the parameters po and p c 
way that p > for any r > 2M. 



(29) 



in such a 
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Figure 1: Accreted fluid velocity u in the linear model (21) versus radial coordinate x. 

On the other hand, model (21) describes the quintessence (not the phantom energy) for the 
entire r range only if p < 0. Consequently, a physically proper description of the quintessence 
can be obtained from (27) if 



For some of the specific choices of a (more specifically, for a = 1/3, 1/2, 2/3, and 1), p(x) 
and u(x) can be calculated analytically (see the Appendix for details on these calculations). In 
Figs. 1 and 2, the radial 4-velocity component, u, and the density normalized to the density at 
infinity, p/poo are plotted against the coordinate x = r/2M. 

Figure 1 (left) shows the plots of u(x) for hydrodynamically stable fluids with c 2 > at 
a = 1/3, 1/2, 2/3, and 1 (the curves are arranged from top to bottom, respectively). Figure 
1 (right) shows the plots of u(r) for c 2 < at a = —1.1, —2, and —1/2 (the curves are also 
arranged from top to bottom). For this case, we chose the boundary condition mh = 1 on 
the horizon. Figure 2 (left) shows the plots of p/poo for a hydrodynamically stable fluid with 
a = 1 for various cases: p = (the model of neutron star matter); po/poo = 16/9 (the linear 
model of nonphantom dark energy); po/Poo — 7/3 (the linear model of phantom energy); and 
Pol Poo = 7/3 (the linear model of phantom energy with pu = 0) (the curves are arranged 
from top to bottom, respectively). Figure 2 (right) shows the plots of (p/poo) for various cases: 
a = —2, po = 0, and A = 4 (the linear model of phantom energy, the upper curve); and 
a = —1/2, po = 0, and A = 4 (the linear model of nonphantom energy, the lower curve). For 
this case, we chose the velocity |uh| = 1 on the horizon. 

3.2 Accretion onto a Moving and Rotating Black Hole 

Let us consider the accretion onto a moving and rotating black hole in the special case of a 
linear equation of state with a = 1. The condition a = 1 allows an exact analytical expression 
to be derived for the accre-tion rate of dark energy onto a black hole. 
For a = 1, we easily find from (23) that 




(30) 



n 



Pcff 



1/2 



n 



Peff.oo 



(31) 
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Figure 2: Accreted fluid density normalized to the density at infinity, p/poc, versus radial 
coordinate x for the linear model (21) (solid curves). The dashed lines indicate the density of 
the A-term, p/poc- 



We obtain the following continuity equation for the particle concentration from (5): 

We can introduce the scalar field <f> in terms of which the fluid velocity can be expressed as 
follows (there is no torsion in the fluid): 

P + P , ( o \ 

n 

We derive an equation for the auxiliary function </> by using Eqs. (31) and (32), 

^ = 0. (33) 

Exactly the same equation arises in the problem of the accretion of a fluid with the equation of 
state p — p [41]. Thus, we reduced the problem of a black hole moving in dark energy with the 
equation of state p = p— po to the problem of a fluid with an extremely hard equation of state, 
p = p. Using the method suggested in [41], we obtain the mass evolution law for a moving and 
rotating black hole immersed in dark energy with the equation of state p = p — p : 

M = An(r 2 + + a 2 )[ Poo + p(poo)}u° Bli , (34) 

where 

r+ = M + (M 2 - a 2 ) 1 ' 2 

is the radius of the event horizon for a rotating black hole, a = J/M is the specific angular 
momentum of the black hole (rotation parameter) , and u ^ is the zeroth 4- velocity component 
of the black hole relative to the fluid. Expression (34) for Ug H = reduces to (18) for a 
Schwarzschild (a = 0) black hole at rest. 



3.3 C haply gin Gas 

Let us consider a Chaplygin gas with the following equation of state as another example of the 
solvable model: 

p = --, (35) 
P 
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where a > 0. The range of parameters p 2 < a represents the phantom energy with a supcrlu- 
minal speed of sound, implying that the phantom energy flux onto the black hole is not fixed 
by the condition of its passage through the critical point. The case of p 2 > a corresponds to 
dark energy with p + p > and < c 2 < 1. We can easily find from Eq. (8) that 



n 



a 



1/2 



The density at the critical point can be calculated from (17) and (36): 

P* = 4 PiL - 3a - 

The velocity and the radial coordinate at the critical point are given by 



ul = 



a 



We then find the constant A from Eq. (9): 

A = A 



a 



4p? 



3/2 



(36) 



(37) 



(38) 



(39) 



For < c 2 < 1, the constant A cannot be smaller than 4, as in the case of the linear model. 
The evolution time scale of the black hole mass without any cosmological change in the dark 
energy density is given by 

(40) 



P 



Note that Eqs. (36)- (39) are applicable only for dark energy with p + p > and are invalid for 
phantom energy On the black hole horizon, 



Ph = jPoo, "H = 



A 
T 



,1/2 



(A/Afpl a 



(41) 



For < c 2 < 1 , the density on the horizon pa cannot be lower than p^ , and «h changes from 1 
to 1/2. The Chaplygin gas density distribution can be determined from the general equations 
(9) and (10): 



3/2 



Q 
72~ 



-1/2 



x A V P 



2 

1 

x 



1/2 



(42) 



The velocity distribution u(r) can be calculated by using Eqs. (9), (36), and (42). In Figs. 3 
and 4, the velocity u and the density normalized to the density at infinity, p/poo, are plotted 
against the coordinate x — r/2M. Figure 3 (left) shows the plots of u(r) for nonphantom dark 
energy at p^/a = 3, 2, and 1.1 (the curves are arranged from top to bottom, respectively). 
Figure 3 (right) shows the plots of u(r) for phantom energy at p^/ot = 0.3, 0.5, and 0.9 (the 
curves are also arranged from top to bottom, respectively). In this case, the boundary con- 
dition «h = 1 is set on the horizon. Figure 4 (left) shows the plots of the normalized density, 
p/poc, for nonphantom dark energy at p^/ot — 3, 2, and 1.1 (the curves are arranged from top 
to bottom, respectively). Figure 4 (right) shows the plots of the normalized density, (p/ Poo), for 
phantom energy at pi^/a — 0.9, 0.5, and 0.3 (the curves are also arranged from top to bottom, 
respectively). For this case, we chose the boundary condition mh = 1 on the horizon. 
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Figure 3: Velocity u versus coordinate x for Chaplygin gas [35] 



4 THE FATE OF BLACK HOLES DURING THE BIG 
RIP 

Let us now consider the evolution of black holes in the cosmological Big Rip scenario, where 
the scale factor a(t) increases to infinity in a finite time [16, 17]. For simplicity, we take into 
account only the dark energy and disregard the other forms of energy. In the linear model (21), 
the Big Rip takes place at p + p < and a < — 1. The following relation can be derived from 
Friedmann's equations in the case of a linear equation of state: 

\p + p\ «a- 3 < 1+a ». 

Setting, for simplicity, p — 0, we find the evolution law of the phantom energy density in this 
Universe: 



t 



where 



3(1 + a) /8tt 



g Poo,i 



1/2 



(43) 



(44) 



Here, poo.i is the initial cosmological phantom energy, and the initial time was chosen in such 
a way that the Big Rip occurs at time r. We easily see from Eqs. (20) and (43) that the Big 
Rip takes place at a = dp /dp < — 1. In general, the condition p + p < alone is not enough 
for the cosmological evolution to be ended with the Big Rip [18]. 

Using Eq. (43), we find the evolution of the black hole mass in the cosmological Big Rip 
scenario from Eq. (18): 

M = Mi (l + —) , (45) 



where 



M t r-t 



M = (3/2)A- 1 |l + a|, 



(46) 



and Mi is the initial mass of the black hole. At a = —2 and a typical value of A = 4 
(which correspond to un = —1), Mq = 3/8. In the limit t — > t (i.e., near the Big Rip), the 
t dependence of the black hole mass becomes linear, M ~ Mq (t — t). When t approaches r, 
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Figure 4: Density normalized to the density at infinity, p/poo, versus coordinate x for model 
(35) (solid curves). The dashed line indicates the normalized density of the A-term, p/ poo- 

the rate of decrease in the black hole mass ceases to depend on the initial black hole mass and 
the phantom energy density: In other words, the masses of all black holes near the Big Rip 
are approximately equal and approach zero. This implies that the accretion of phantom energy 
dominates over the Hawking evaporation until the black hole mass decreases to the Planck mass. 
Formally, however, all black holes in the Universe completely evapo-rate during the Hawking 
radiation in the Planck time before the Big Rip occurs. 

5 THE ACCRETION OF A SCALAR FIELD 

In this section, we compare our calculations of the accretion of an ideal fluid with similar 
calculations of the accretion of a scalar (nonphantom) field onto a black hole [24-29]. The dark 
energy is commonly modeled by a scalar field with a potential V{<f>). The approximation of an 
ideal fluid is rougher, since the scalar field <f) and d^cfi cannot be unambiguously reproduced for 
given p and p, which characterize an ideal fluid. Despite this difference between the scalar field 
and the ideal fluid, we will show that our results are in close agreement with the corresponding 
calculations of the accretion of a scalar field onto a black hole. 

The Lagrangian of the scalar field is L = K — V , where K is the kinetic term and V is the 
potential. For the standard choice of the kinetic term 

K = cj>,^/2, 

the corresponding energy flux onto the black hole is 

Tor = 4>,t4>,r- 

Jacobson [24] found a solution for the scalar field in the Schwarzschild metric for a zero potential, 
V = 0: 

= ^ oo [t + 2Mln(l-2M/r)], 

where (poo is the scalar field at infinity. Frolov and Kofman [26] showed that this solution is 
also valid for many scalar fields with a nonzero potential V() under certain conditions. For this 
solution 

TJ = -(2M) 2 ^/r 2 , 
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and, accordingly, 

M = 4tt(2M) 2 ^. 

The energy-momentum tensor constructed using Jacobson's solution is identical to the 
energy-momentum tensor for an ideal fluid with an extremely hard equation of state, p = p, 
after the substitution 

Poo -» </>L A Poo -» ^L/ 2 - 

This is not surprising, since the theory of a scalar field with a zero potential, V(<f>), is identical 
to the model of an ideal fluid [48]. In view of this correspondence, we easily see agreement 
between our result (18) for M in the case of p = p and the corresponding results from [24, 26]. 

The Lagrangian of the scalar field that describes the phantom energy must have a negative 
kinetic term [16, 17], for example, 

(see [49] for more general cases) . In this case, the phantom energy flux onto the black hole has 
the opposite sign, 

To r = —<t>,t<t>,r, 

where 4> is the solution of the same Klein-CGordon equation as that for the standard scalar 
field, but with the substitution V — > — V. For a zero potential, this solution is identical to 
Jacobson's solution [24] obtained for a scalar field with a positive kinetic term. 

However, the Lagrangian with a negative kinetic term and V(<j>) = does not describe the 
phantom energy. At the same time, the solution for a scalar field with V(<f>) = is identical to 
the solution for a positive constant potential, V = const, which can be chosen in such a way 
that 

p = -j> 2 /2 + V > 0. 

In this case, the scalar field describes the accreted phantom energy with p > and p < —p, 
which leads to a decrease in the black hole mass at the rate 

M = -M2M) 2 4>lo- 

A simple example of phantom cosmology (but without the Big Rip) is realized by a scalar 
field with the potential 

V = toVA 

where m ~ 10 -33 eV [50]. After a short transition period, this cosmological model approaches 
an asymptotic state with 

H ~ m0/3 1/2 , </»~2m/3 1/2 . 

In the Klein-Gordon equation (with the substitution V — > — V mentioned above), the term m 2 
becomes equal to the other terms only on the scale of the cosmological horizon, implying that, 
in this case, Jacobson's solution is also valid. Calculations of the corresponding energy flux 
onto the black hole yield 

M = -4tt(2M) 2 2 c = -64M 2 m 2 /3. 
For Mo = Mq and m = 10~ 33 cV, the effective time of the decrease in black hole mass is 

r = (3/64)M- 1 m- 2 - 10 32 yr. 
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6 DISCUSSION AND CONCLUSIONS 

In recent years, the concept of dark energy has been accepted and extensively discussed in 
cosmology The possible existence of dark energy with a negative pressure leads to new cosmo- 
logical scenarios, including the exotic model of the Universe in which all of the bound objects 
are destroyed and which dies itself as a result of the Big Rip. To determine the fate of black 
holes in this cosmological scenario, we considered the spherically symmetric, stationary accre- 
tion of dark energy mod-eled by an ideal fluid onto a black hole. We derived general equations 
for the accretion of an ideal fluid with the equation of state p = p(p) onto a Schwarzschild black 
hole. In particular, these equations can be used to describe the accretion of thermal radiation, 
dark energy, and phantom energy. We also considered the accretion onto a moving and rotating 
black hole in the special case of an extremely hard equation of state, p = p. We calculated the 
change in the black hole mass through accretion. The black hole masses for p+p > were found 
to increase, as in the usual case. However, a qualitatively new result was obtained for phantom 
energy, i.e., for a medium with p + p < 0. We found that the black hole masses decrease in 
this situation. Using this result, we solved the problem of the fate of black holes in a universe 
that undergoes the Big Rip. It turns out that all black holes in this Universe must decrease 
their masses and disappear completely by the Big Rip. We also considered the correspondence 
between the accretion of dark energy in the model of an ideal fluid and the accretion of a scalar 
field. 

Acknowledgments: This work was supported in part by the Russian Foun-dation for Basic 
Research, project nos. 02-02-16762a, 03-02-16436a, and 04-02-16757a and by the Ministry of 
Science of the Russian Federation, grants 1782.2003.2 and 2063.2003.2. 



Appendix: Analytical Solutions for p(x) and u(x) 

In the model under consideration (Section 3), analytical solutions can be found for the de- 
pendence of the dark energy density and accretion rate on radius r. Using Eq. (23) for the 
concentration and Eq. (11) for the constant C2, we derive the following equation for p c g from 
Eqs. (9) and (10) 



PcS 
Peff.oo 



2a/(l+a) 



2 A 2 



X X \/?cff,oo 



PcS 



-2/(l+a)' 



1. 



Defining 



y 



PcS 



,Peff, 00 

we obtain the following equation from (47): 



2/(l+a) 



y 1-- 



y 



l-a 



A 2 

-IT =0, 



(47) 



(48) 



(49) 



which can be solved analytically for certain values of a. For a = 1/3, Eq. (49) reduces to a 
cubic equation: 

2\ o A 2 



1 - 



z 2 + 



X J XT 

where z = y 1 ^ 3 . Solving this equation yields the fluid density distribution for a = 1/3: 

i 2 



P = 



Po 



(*•-?) 



z + 



1 

3(1 - 2a;" 1 ) 



(50) 



(51) 



13 



where 

>,/f cos ( 4f - f V 2 < a; < 3, 

z={ >A ; - - ' (52) 

2^/f cos (§ J , a: > 3, 



/3 = arccos 

and 



2 (a/3) 3 / 2 



(53) 

108 



a = 9i & = T - 7^ ~ ; x 7 . (54) 

3(1 -2/x) 2 27(1 -2/x) 3 (l-2/x)x 4 

This solution corresponds to a thermalized photon gas in which the photon mean free path 
is much smaller than the radius of the black hole horizon, A/ p <C 2M. In this situation, the 
photon gas may be treated as an ideal fluid. In the opposite case, Xf p ^> 2M, the photons are 
free particles, and their accretion rate is determined by the well-known cross section for the 
gravitational cap-ture of relativistic particles by a black hole. The corresponding accretion rate 
is 

M = 27ttM 2 Poo 

The case of a = 2/3 is similar to the case considered above. We obtain the following 
equation instead of (50): 

*>(l-l)- z+ *=0, (55) 



x , 

where again z — y 1 / 3 . The fluid density distribution in this case is 



/' T/V • ( Poo ilh, ) z 5/2 , ( •">() 



where z is given by 



[3 is defined by Eq. (53), and 



cos (?f - f ) , 2 < x < 9/4, 
2yfcos(§), x>9/4, 

1 , 2187 V3 



(57) 



1-2/x' 128 (1 - 2/x) x 4 ' 

For a = 1/2, (49) is a quadratic equation and has a simple analytical solution: 

p = f + ("» - ?) 23/2 ' (58) 



where 

[1-3125 (1-2/x) (16a; 4 )- 1 ] 172 } (1-2/x)" 1 , 2 < x < 5/2 
[1-3125 (l-2/x)(16x 4 

For a = 1, Eq. (49) is linear in y, which gives 



±<jl + [1-3125 (1-2/x) (16X 4 )- 1 ] 172 } (1-2/x)" 1 , x > 5/2. 
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